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Abstract 

I discuss the theoretical uncertainties in the extraction of | V c b\ from 
semileptonic decays of B mesons, taking into account the most recent the- 
oretical developments. The main sources of uncertainty are identified both 
for the exclusive decay mode B — > D*tv and for the inclusive channel 
B —* Xlv. Prom an analysis of the available experimental data, I ob- 
tain | V^lcxd = 0.041 ± 0.003 CX p ± 0.002 t h from the exclusive mode, and 
I Kfelinci = 0.040 ± O.OOlexp ± 0.005 t h from the inclusive mode. I also give 
a prediction for the slope of the form factor T{yS) at zero recoil, which is 
q 2 = 0.8 ±0.3. 
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1 Introduction 



Semileptonic decays of B mesons have received a lot of attention in recent years. 
The decay channel B — > D*£ v has the largest branching fraction of all B-meson 
decay modes, and large data samples have been collected by various experimental 
groups. From the theoretical point of view, semileptonic decays are simple enough 
to allow for a reliable, quantitative description. Yet, the analysis of these decays 
provides much information about the strong forces that bind the quarks and 
gluons into hadrons. Schematically, a semileptonic decay process is shown in 
Fig. |I[ The strength of the b — > c transition vertex is governed by the element 
V c b of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The parameters of this 
matrix are fundamental parameters of the Standard Model. A primary goal of 
the study of semileptonic decays of B mesons is to extract with high precision 
the values of V c b and V u b- The problem is that the Standard Model Lagrangian 
describing these transitions is formulated in terms of quark and gluon fields, 
whereas the physical hadrons are bound states of these degrees of freedom. Hence, 
an understanding of the transition from the quark to the hadron world is necessary 
before the fundamental parameters can be extracted from experimental data. 




q 

Figure 1: Semileptonic decay of a B meson. 

Today, our knowledge of the elements of the CKM matrix, as extracted 
from direct measurements of flavour-changing transitions, is as follows: The best 
known entries are V u d and V us , which have an uncertainty of 0.1% and 1%, re- 
spectively. The next well known entry is already related to the 6-quark; V c b is 
now determined to an accuracy of 7%. Then follow V c d, V cs , and V u b, with 10%, 
20%, and 30% uncertainties, respectively. No direct measurements exist for the 
matrix elements related to the top quark. 

In this talk I discuss the status of the theoretical developments underlying the 
determination of V c b, both from exclusive and from inclusive semileptonic decays 
of B mesons. 
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2 \Vcb\ from Exclusive Decays 

With the discovery of heavy-quark symmetry (for a review see Ref. and ref- 
erences therein), it has become clear that the study of the exclusive semileptonic 
decay mode B — > D*iv allows for a reliable determination of | V c b\, which is free, 
to a large extent, of hadronic uncertainties @— 1|. Model dependence enters this 
analysis only at the level of power corrections of order (Aqcd/^q) 2 -[| These cor- 
rections can be investigated in a systematic way, using the heavy-quark effective 
theory |J. They are found to be small, of order a few per cent. 

The analysis consists in measuring the recoil spectrum in the decay B — > 
D*£ v. One introduces the kinematic variable 

E D * m% + m%, - q 2 

w = v B ■ v D * = = , (1) 

m D * 2m B m D * 

which is the product of the four-velocities of the mesons. Here Ed* denotes the 
recoil energy of the D* meson in the parent rest frame, and q 2 = (ps — Pd*) 2 is 
the invariant momentum transfer. The differential decay rate is given by {|, ||| 

' (m B — run*) 2 m£)*Vw 2 — 1 (w + l) 2 



dw 48ir 3 

Aon ™2 _ „„™_™_ _l™2 1 

V cb \ 2 F 2 (w). (2) 



4w m 2 B — 2wmBiriD* + m 2 D , 



w + 1 (rriB — m D *) 2 

The function ^(w) denotes the (suitably defined) hadronic form factor for this 
decay. It is conventional to factorize it in the form J-'iw) = t}a £(w), where t]a 
is a short-distance coefficient, and the function £(w) contains the long-distance 
hadronic dynamics. Apart from corrections of order Aqcd/^q, this function 
coincides with the Isgur-Wise form factor |3], |7J . In analogy to the case of light- 
quark SU(3) flavour symmetry, in which the Ademollo-Gatto theorem protects 
the K — > 77 transition form factor against first-order symmetry-breaking correc- 
tions at q 2 = P|, there is a theorem which protects the function £(w) against 
first-order Aqcd/^q corrections at the kinematic point of zero recoil (w = 1). 
This is Luke's theorem 0, which determines the normalization of £(u>) at w — 1 
up to corrections of order (Aqcd/^q) 2 , i-e. ^(1) = 1 + 6y m 2. 

The strategy for a precise determination of | V c b\ is thus to extract the product 
| V c b \ ^F(w) from a measurement of the differential decay rate, and to extrapolate 
it to w — 1 to measure 

| Vcb\ .F(l) = | Vcb\ T) A (1 + Si/ m a) . (3) 

The task of theorists is to provide a reliable calculation of the quantities T]a and 
8u m 2 in order to turn this measurement into a precise determination of | V c b\- I 
will now discuss the status of these calculations. 



1 I shall use Aqcd ~ 0.25 GeV as a characteristic low-energy scale of the strong interactions, 
and ttiq as a generic notation for m c or TOf,. 
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2.1 Perturbative corrections 



The short-distance coefficient i]a takes into account the finite renormalization of 
the axial vector current arising from virtual gluon exchange. It can be calculated 
in perturbation theory. At the one-loop order, one finds JTU], 0, [TTJ 

a s (M) fm b + m c m b 8\ 
VA = 1 + In . 4 

7r \rrii, — m c m c 3/ 

The scale M of the running coupling constant is not determined at this order. 
A reasonable choice is to take the average virtuality of the gluon in the one-loop 



diagrams In the case of t]a, this so-called BLM scale has been calculated 
to be M = 0.51y / mi/n\. (in the MS scheme) |R|. Taking then mj = 4.8 GeV, 
fn c / m b = 0.30 ± 0.05, and Aq C d = 0.25 GeV for the scale parameter in the 
two-loop expression for the running coupling constant, one obtains values in the 
range r] A = 0.950-0.965. 

Several (partial) higher-order calculations have been performed to improve 
this result. Using renormalization-group techniques, logarithms of the type (a s In z) r 
a s (a s In z) n , and (m c /mb)(a s In z) n , where z = m c /mb, have been resummed to 
all orders in perturbation theory [|T4 — [ 18 1 . This leads to the somewhat larger 



value t]a — 0.985. Another class of higher-order corrections consists of the so- 
called renormalon chain contributions, which are terms of order I3q~ 1 o.™ in the 
perturbative series for t]a- Resumming these terms to all orders gives the lower 



value r] A ~ 0.945 19 



The main virtue of these partial higher-order calculations is to provide an 
estimate of the theoretical uncertainty in the value of r] A . Thus, I quote the final 
result as 

T] A = 0.965 ±0.020. (5) 



2.2 Power corrections 

Hadronic uncertainties enter the determination of | at the level of second- 
order power corrections, which are expected to be of order (AqcdA^c) 2 ~ 3%. 
For a precision measurement, it is important to understand the structure of these 
corrections in detail. This is the most complicated aspect of the theoretical anal- 
ysis, which unavoidably introduces some amount of model dependence. However, 
since the goal is to estimate an effect which by itself is very small, even a large 
relative error in 5 x / m 2 is acceptable. 

Three approaches have been suggested to estimate these corrections. The idea 
of the "exclusive" approach of Falk and myself is to classify all 1 / mg operators in 
the heavy-quark effective theory and to estimate their matrix elements between 
meson states p0| . This last step is model-dependent. A typical result obtained 
in this way is 8x/ m 2 = —(3 ± 2)%. In Ref. [|IJ, the error has been increased to 
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±4% in order to account for the model dependence and unknown higher-order 
corrections. A similar result, —5% < 5i/ m a < 0, has been obtained by Mannel 







The idea of the "inclusive" approach of Shifman et al. is to apply the operator 
product expansion to the 5-meson matrix element of the time-ordered product of 
two flavour-changing currents, and to equate the resulting theoretical expression 
to a phenomenological expression obtained by saturating the matrix element with 



physical intermediate states [22]. This leads to sum rules that imply inequalities 
for the B — > D* transition form factors. In particular, one obtains the bound 
fii/m 2 < ~ | ~ m B)/ m c — — 3%. The authors of Ref. |2^[ make an "educated 
guess" that the value of Sy m 2 is actually much larger, namely — (9 ± 3)%. 



It is possible to combine the above predictions in a "hybrid" approach, which 
uses the sum rules to put bounds on the hadronic parameters that enter the "ex- 
clusive" analysis || . One then finds that for all reasonable choices of parameters 
the results are in the range —8% < 5 ly / m 2 < —3%, which is consistent with all 
previous estimates at the lex level. Thus, I quote the final result as 

5 1/m2 = -(5.5 ±2.5)%. (6) 



2.3 Determination of \ V c b\ 

Combining the above results, I obtain for the normalization of the hadronic form 
factor at zero recoil 



/m 2 



0.91 ±0.04. 



(7) 



To be conservative, I have added the theoretical errors linearly. Three exper- 
iments have recently presented new measurements of the product |Vd,|.F(l). 
When rescaled using the new lifetime values t b o = (1.61 ± 0.08) ps and t b + = 
(1.65 ±0.07) ps 



, the results are 



0.0351 ±0.0019 ±0.0020 
V cb \ = <j 0.0364 ± 0.0042 ± 0.0031 
0.0388 ±0.0043 ±0.0025 



Ref. [24], 
Ref. IE25I, 



(8) 



Ref. 26 



where the first error is statistical and the second one systematic. I will follow 
the suggestion of Ref. []27| and add 0.001 ± 0.001 to these values to account for 
a small positive curvature?] of the function ^F(w) . Taking the weighted average 
of the experimental results, which is | V c b\ = 0.0370 ± 0.0025, and using the 
theoretical prediction (|7|), I then obtain 



cb\ 



0.0407 ± 0.0027 exp ± 0.0016 th 



(9) 



2 The fit results are obtained assuming a linear form of F{i 
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which corresponds to a measurement of \V c b\ with 7% accuracy. 
2.4 Prediction for the slope parameter g 2 

In the extrapolation of the differential decay rate (0) to zero recoil, the slope 
of the function J-(w) close to w = 1 plays an important role. One defines a 
parameter g 2 by 

f(w) = F(l){l-g 2 (w-l) + ...}. (10) 

It is important to distinguish g 2 from the slope parameter q 2 of the Isgur- 
Wise function. They differ by corrections that break the heavy-quark symmetry. 
Whereas the slope of the Isgur-Wise function is a universal, mass-independent 
parameter, the slope of the physical form factor depends on logarithms and in- 
verse powers of the heavy-quark masses. The relation between the two parameters 
is§ 

f = g 2 + (0.16 ±0.02) + 0{l/m Q ). (11) 

An estimate of the Aqcd / rnq corrections to this relation is model-dependent and 
thus has a large theoretical uncertainty. I shall not attempt it here. 

The slope parameter of the Isgur-Wise function, g 2 , is constrained by the 
Bjorken [p8| , |29| and Voloshin |Kj sum rules. At the tree level, it was known for a 



long time that 1/4 < g 2 < pa 0.75. However, only recently Grozin and Korchem- 
sky have shown how to include perturbative and nonperturbative corrections to 
these bounds |32|. The results are shown in Fig. 0. Here, the scale parameter 



\i has to be chosen large enough for the operator product expansion to be well 
defined, but it is otherwise arbitrary. Assuming that values fi > 0.8 GeV are 
sufficiently large, one finds that g 2 is constrained to be very close to 0.6. This 
value is in good agreement with earlier predictions obtained from QCD sum rules, 
which gave g 2 = 0.7 ± 0.1 |, 0, 0. 

From Fig. |[ and using ([□]), I conclude that 

g 2 = 0.65 ±0.15, g 2 = 0.8 ±0.3. (12) 

This prediction compares well with the average value observed experimentally, 
which is g 2 = 0.87 ± 0.12 [El-pJ. 



3 \Vcb\ from Inclusive Decays 

Complementary to the analysis of exclusive decays is the study of the inclusive 
semileptonic decay rate for B — > X tv. Since | V u b/V c b\ 2 < 1%, one can to very 
good approximation neglect the contribution of charmless final states and consider 
X to be a hadronic state containing a charm particle. An obvious advantage 
of inclusive decays is the existence of high-statistics data samples. From the 
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Figure 2: Bounds for the slope parameter q 2 following from the Bjorken 
and Voloshin sum rules, from Ref. 



32] . The point with the error bar 



shows the QCD sum rule prediction. 



theoretical point of view, summing over many final states eliminates part of the 
hadronic uncertainty. 

As in the exclusive case, the framework for the theoretical description of 
inclusive decays is provided by the heavy-quark expansion. It could be shown that 
the leading term in this expansion reproduces the free-quark decay model, while 
the nonperturbative corrections to this model can be systematically included in 
an expansion in powers of l/m b [35|-|41|. The total semileptonic decay rate can 
be written as 



T = -^L\ Vcb \ 2 ml 

192vr 3 I cfc l b 



1 + 



Xi + 3A 2 
2m\ 



f(m c /m b ) 



6A, 



m 



\ g(m c /m b ) + 



71 



mi 



(13) 



where the ellipsis represents terms of higher order in l/m b or a s . In this expres- 
sion, m b and m c denote the pole masses (defined to the appropriate order in per- 
turbation theory) of the heavy quarks, f(m c /m b ) ~ 0.52 and g(m c /m b ) ~ —0.87 
are kinematic functions, and Ai and A2 are nonperturbative hadronic parameters. 
I will now discuss the theoretical uncertainties in the evaluation of (O). 



3.1 Perturbative corrections 

Let me first discuss the uncertainty due to unknown higher-order perturbative 
corrections. Only the correction of order a s is known exactly P2| . However, 
recently Luke et al. have computed the part of the next-order term that depends 
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on the number of light-quark flavours | |4"3|| . In the MS scheme, the result is 

r = t _ x 67 ^K)_ (1 68 A+> /^K)A 2 + _ = i_o.n_o.Q6_. . . (14) 

where To is the decay rate at the tree level, and Po = 11 — |rij is the first 
coefficient of the /5-function. If one uses this partial calculation to estimate the 
uncertainty, which for an asymptotic series is given by the size of the last term 
to be included, one finds (ST/T) pert ~ 6%, which is in good agreement with an 
estimate of the renormalization-scale and -scheme dependence by Ball and Nierste 
M\. Recently, Ball et al. have performed an all-order resummation of the terms 



of order (5q ~ 1 a™ for the above series ||45| . They find that the effect of higher-order 
terms is important and leads to r/T = 0.77 ± 0.05. Note that this corresponds 
to an effective scale M ~ 1 GeV in (|l"3"D, which is rather low. The corresponding 
value of the coupling constant is a s (M) ~ 0.43. It is difficult to derive a reliable 
error estimate from these analyses, but I think a reasonable number is 



5T 



) -10%. (15) 



pert 



3.2 Power corrections 

The leading power corrections in the expression for the inclusive decay rate appear 
at order 1/m^. They are proportional to two hadronic parameters with a simple 
physical interpretation: Ai is related to the average momentum of the 6-quark 
inside a B meson at rest, and A2 is proportional to the vector-pseudoscalar mass 
splitting. I shall use 



2 



\ 1 = -(p b 2 ) = -(0.4 ±0.2) GeV 
m|*-m| ^ Q12 Gey2 
4 

The value of Ai is a compromise between the theoretical estimates obtained in 



Refs. [46, 47 



The power corrections reduce the total decay rate by — (4.2±0.5)%, which is a 
rather small effect. The uncertainty in the value of Ai introduces an uncertainty 
on the value of T of 0.6%, which is almost negligible. I increase this value in 
order to account for higher-order power corrections, which I expect to be of order 
(Aqcd/^c) 3 ~ 0.5%, and quote 

^r) ^1%. (17) 

J- / power 

This is the smallest contribution to the total theoretical uncertainty. 
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3.3 Dependence on quark masses 



Another source of nonperturbative uncertainty results from the appearance of 
the heavy-quark masses in the expression for the inclusive decay rate. The pole 
masses of the bottom and charm quarks have an uncertainty of at least several 
hundred MeV. Since the rate is proportional to m|, this seems to be a severe 
limitation. However, it has been pointed out that the actual uncertainty is lower 
due to a strong correlation between the values of the two heavy-quark masses 



2^| . In fact, one should consider the decay rate as a function of nib and of 
the difference Am = nib — m c . I shall assume that nib has an uncertainty of 
300 MeV. However, the mass difference is known to much higher precision. Using 



the heavy-quark expansion, one can derive that p0 | 

Am = m b — m c = (m B — m D ) {l — ^ — h 0(llrn? o ) 

I 2m B m D v 



where ijib = \ (mg + 3ms*) = 5.31 GeV and mo = \ + 3ms*) = 1.97 GeV 
denote the spin-averaged meson masses. This relation leads to 

Am = (3.40 ± 0.03 ± 0.03) GeV , (19) 

where the first error reflects the uncertainty in the value of Ai, and the second one 
takes into account unknown higher-order corrections. Hereafter, I shall assume 
an uncertainty of 60 MeV in the value of Am. 

In Fig. |3|, I show the dependence of the decay rate on these two param- 
eters, using the value a s (M) = 0.4 for the strong coupling constant in (JT3|). 
Clearly, the variation with Am is much stronger than the variation with mj. 
For rrib = 4.8 GeV and Am = 3.4 GeV, I find the partial derivatives <5r/r ~ 
5.73 5(Am)/Am and 5T /T ~ — 0.55 5mb/nib. Since the errors in m;, and Am are 
essentially uncorrected, this leads to 



r y massfiS V v 60 Mevy v 300 Mev; v ; 



Note that this is dominated by the (rather small) uncertainty in the mass differ- 
ence Am. 



3.4 Determination of \V c b\ 

Adding the above errors linearly and taking the square root, I conclude that the 
theoretical uncertainty in the extraction of | V c b\ from inclusive decays is 

T^T-11%. (21) 

I Vcb\ 
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From an analysis of the experimental data, one then obtains |27 



■ i 0.039 ± 0.001 cxp ± 0.005th ! measurements at T(4s), , 
1 ' /,! = I 0.042 ± 0.002 exp ± 0.005 th ; measurements at Z°. 1 ' 



The theoretical uncertainty in these numbers is somewhat larger than in the case 
of the exclusive analysis; however, the experimental errors are smaller. 



4 Summary 

The most precise measurements of the element V c b of the CKM matrix come from 
the analysis of semileptonic decays of B mesons. From the measurement of the 
recoil spectrum in the exclusive channel B — > D*£ u, one obtains 

I Kilexd = 0.041 ± 0.003 cxp ± 0.002 th , (23) 

where the theoretical error is dominated by the uncertainty in the calculation of 
nonperturbative power corrections of order (Aqcd/^q) 2 - From the measurement 
of the total inclusive decay rate, on the other hand, one finds 

I V^linci = 0.040 ± O.OOlexp ± 0.005 th . (24) 

In this case, the main theoretical uncertainty comes from the uncertainty in 
the value of the quark mass difference — m c , as well as from higher-order 
perturbative corrections. 



9 



Given the fact that both methods are very different both from the experimen- 
tal and from the theoretical point of view, it is most satisfying that the results 
are in perfect agreement. Combining them, I obtain the final value 

\V cb \ = 0.041 ±0.003. (25) 
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